Length-Controlled Elasticity in 3D Fiber Networks 
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We present a model for disordered 3D fiber networks to study their linear and nonlinear elasticity 
over a wide range of network densities and fiber lengths. In contrast to previous 2D models, these 
3D networks with binary cross-links are under-constrained with respect to fiber stretching elasticity, 
suggesting that bending may dominate their response. We find that such networks exhibit a fiber 
length-controlled bending regime and a crossover to a stretch-dominated regime for lengths beyond 
a characteristic scale that depends on the fiber's elastic properties. Finally, by extending the model 
to the nonlinear regime, we show that these networks become intrinsically nonlinear with a vanishing 
linear response regime in the limit of floppy or long filaments. 
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Materials ranging from paper and textiles to the struc- 
tural components of living cells and tissues [T] all exhibit 
networks of fibers or stiff polymers. Such networks have 
extraordinary mechanical properties [SHI]. Their elastic- 
ity depends in part on their connectivity [FJ |Bj, in analogy 
with jammed matter[7l[8] and random network glasses[9]- 
The mechanics of the constituent fibers, and specifi- 
cally their bending rigidity can also strongly impact net- 
work elasticity. However, the relative importance of fiber 
stretching versus bending is not understood, especially in 
3D. Prior work has mostly focused on 2D networks 1 1 QffT5] 
since simulations in 3D have been proven to be challenging 
and have usually been limited to small system size|16l 
ITS] , Significant qualitative differences are expected be- 
tween 2D and 3D networks: for the typical case of binary 
fiber interactions, the high- molecular weight limit in 2D 
actually corresponds to the Maxwell central-force (CF) 
isostatic threshold, where stretching interactions begin to 
completely constrain network deformations. In contrast, 
3D networks with binary interactions remain well below 
CF isostaticity. Thus, owing to their marginal stability, 
real 3D fiber networks are expected to be fundamentally 
more bending-dominated and more prone to collective 
nonaffine deformations |16l 118] , 

Here we develop a numerical model for the elasticity of 
random 3D fiber networks with binary cross-links. This 
model provides access to a wide range of network densi- 
ties below the CF isostatic threshold, as well as the previ- 
ously inaccessible high-molecular weight limit. These net- 
works exhibit various qualitatively distinct linear elastic 
regimes: a critical regime governed by the rigidity perco- 
lation point, a length-controlled bending regime and an 
affine stretching regime, as illustrated in Fig. [lji. We pro- 
vide a scaling analysis for insight into the origins of these 
regimes. Network mechanics is determined by the ratio 
between a nonaffinity length scale and molecular weight. 
The high-molecular weight limit exhibits behavior rem- 
iniscent of zero-temperature critical behavior, including 
divergent strain fluctuations; these ultimately govern a 
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nonlinear elasticity (b) of 3D fiber networks on the Phantom 
FCC lattice, where L is the average filament length, z is net- 
work connectivity, 7 is strain and « is the fiber bending rigid- 
ity. All lengths are measured in units of the lattice spacing Iq 
and K in units of /i£q. 

crossover from the length-controlled bending regime to an 
affine, stretching dominated network response. This is re- 
markable, since the network remains well below Maxwell's 
CF isostatic connectivity threshold in this limit. Thus, 
paradoxically, although such networks can only be rigid 
at non-zero fiber bending stiffness, we find that no mat- 
ter how weak this bending rigidity is, the network elas- 
ticity becomes insensitive to fiber bending in the limit 
of high — yet finite — molecular weight. Moreover, in the 
limit of floppy filaments with weak bending rigidity or 
high molecular weight, these networks become intrinsi- 
cally nonlinear with a vanishing linear response regime 
(Fig#). 

Much has been learned about stiff polymer gels from 
minimal models, such as 2D Mikado networks of randomly 
placed straight filaments with binary cross-links [W[ ITT]. 
The elasticity of such Mikado networks is governed 
by nonaffine fiber bending deformations at low densi- 
ties, while higher density networks exhibit predominantly 
affine stretching elasticity of single fiber segments. This 
nonaffine- affine (NA-A) transition can be understood as 
being the result of increasing fiber-length. However, for 
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FIG. 2: (Color online) a) The Shear modulus as a function 
of L in units of £o for various k in units of /^q. Co is the 
(affine) shear modulus of the undiluted Ph-FCC lattice. The 
inset illustrates the phantom principle: At each lattice vertex 
3 independent binary cross-links are formed between randomly 
chosen fiber pairs labeled by color, b) Non-affinity parameter 
r as a function of L. Dashed black lines indicate a slope of 2. 

such 2D networks, this high molecular weight limit actu- 
ally coincides with Maxwell's CF isostatic connectivity, 
2 CF = 2d in d dimensions [5], which can also give rise 
to a NA-A transition 6J; it is thus unclear whether the 
observed transition in 2D is controlled by CF stretching 
constraints, or by filament length, as previously suggested 
by scaling arguments and floppy mode theory [TDHT^] . 
However, 3D networks with binary cross-links — like most 
biopolymer systems — are qualitatively different; in this 
case the high-molecular weight limit corresponds to net- 
work connectivities well below z CF . In the absence of fiber 
bending resistance, such networks exhibit zero-energy de- 
formation modes and hence, they do not resist shear 
stresses. Thus, there are reasons to question the existence 
of an affine limit in realistic 3D networks with fibers that 
are softer to bending than to stretching |12l IT51 [TO] . This 
is still subject of debate since studies in 3D have so far 
been limited to small systems [TH] or to networks with high 
connectivities [19] . 

To provide insight in network mechanics at densities 
ranging from the rigidity percolation point to the high- 
molecular weight limit, we develop a 3D lattice-based 
fiber network model with binary cross-links. Our net- 
work's consists of straight fibers organized geometrically 
on a face centered cubic (FCC) lattice. However, we limit 
the maximum coordination number to four by randomly 
choosing three independent pairs of cross-linked fibers out 
of the six fibers crossing at every vertex. Although the dif- 
ferent binary cross-links may overlap geometrically, they 
do not constrain each other [2T] (inset Fig. [2^,). Therefore, 



we term this the Phantom FCC (Ph-FCC) lattice. This 
model is similar to a generalized Kagome lattice in 3D [20] . 
although it has a lower symmetry than the Ph-FCC lat- 
tice. By randomly cutting bonds with a probability, 1— p, 
we tune the average molecular weight, L = £q/(1 — p), 
where £o is the distance between vertices [SI [21], 

The elastic energy of the 3D Ph-FCC network involves 
stretching and bending contributions of the constituent 
fibers, characterized by their stretching modulus // and 
bending rigidity n. Each lattice vertex consists of 3 in- 
dependent freely-hinging binary cross-links ranked by h. 
For small displacements, denoted by uf, the stretching 
energy of the network is expressed as 

h=i (ij) 

where the second sum extends over neighboring pairs of 
vertices, = — u' 1 and f ^ is the bond direction in 
the undcformed lattice. Bond-dilution is implemented by 
setting g^j = 1 for present bonds and g^ = for removed 
bonds. Fibers form straight chains that resist angular 
deflections, leading to a total bending energy [6) [14] . 

£; B = ^EE4^[K-u^)xr lJ ] 2 . (2) 

h=l (ijk) 

Since the cross- links themselves do not provide a torsional 
stiffness, the second sum only extends over coaxial nearest 
neighbor triplets along the same fiber. The shear mod- 
ulus, G, is determined numerically by applying a shear 
strain along the Ill-plane with Lees-Edwards periodic 
boundary conditions and energy minimizations are per- 
formed using a conjugate gradient algorithm. 

Crucially, to avoid effects due to filaments spanning the 
network, thereby making unphysical stretch contributions 
to the elasticity of the sample, at least one bond is re- 
moved along every fiber. Our network sizes range from 
W z = 20 3 to 150 3 unit cells, with up to three times that 
many cross-links. Due to the finite system size, this model 
can only approach z = 4 asymptotically from below. 

Linear regime We find that these networks have a fi- 
nite shear rigidity only if n > 0, even though the perfect 
(z = 4) Ph-FCC lattice deforms affinely and has a fi- 
nite shear modulus for k — 0. Nonetheless, over a broad 
range of L, the system can be bending dominated G ~ K 
(low k), or stretching dominated, G ~ /i (high k), as 
shown in Fig. [2^. Interestingly, there appear to be two 
distinct regimes well above the rigidity percolation point: 
a bending-dominated regime where G depends strongly 
on L (low k and L) and an L- and K-independent affine 
stretching regime (high k and L). 

These results can be understood as follows. In the high- 
k limit, the system deforms affinely, with a shear modulus 
G A ~ jlsZ. However, in the critical regime — controlled 
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by the bending rigidity percolation point Zb — G vanishes 
continuously with Az = z — zt,[6l [HI EH] as 
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for high and low k, respectively. We find Zb ps 2.4 and 
/ ~ 0.65 for a system size VF 3 = 30 3 , as demonstrated 
in the inset in Fig. [3]by showing that G\Az\~t / k reaches 
a plateau for low values of Az. The rigidity threshold is 
similar to observations in prior 3D models[18], although 
/ is considerably lower here, which is more consistent 
with findings on the generalized 3D Kagome lattice [213], 
The rigidity threshold can be estimated using a Maxwell 
counting argument 5, 6, 1 1 8 j : this connectivity-threshold 
occurs when per cross-link the number of stretching con- 
straints, and bending constraints, n^ld — l)z 2 /16, 
equal the number of internal degrees of freedom d. Here, 
the number of bonds per cross-link rib — 2 in the undi- 
luted network (z = 4). This yields Zb rj 2.6, in reasonable 
agreement with the numerical results. 

Since the CF isostatic point lies far beyond the physical 
connectivity range of this model, the naive expectation 
would be that this percolation regime extends over the 
whole range z < 4. This would imply bending dominated 
network elasticity for low k, such that G <C Ga as z — >• 4 
from below, with a discontinuous transition at z = 4 to 
affine stretching dominated behavior due to fibers that 
span the whole network. However, this argument ignores 
possible effects due to filament length. In networks of 
straight fibers with binary interactions, the average fiber 
length diverges as z —¥ 4 and large L may lead to nonaffine 
displacements over greater length scales [12]. The effects 
of high L on the deformation field have been discussed 
in the context of 2D Mikado networks using both scaling 
arguments [TO] ITT] and floppy mode theory [12]. 

Here, we investigate the effects of molecular weight 
on the nonaffine deformation field and their implications 
for the mechanics of 3D fiber networks. Network nodes 
along a fiber can undergo independent nonaffine defor- 
mations scaling as jL to avoid stretching of the other 
fibers to which they are connected. This direct scaling 
of nonaffine displacements with L constitutes one of the 
central assumptions of floppy mode theory that was ap- 
plied to Mikado networks [12]. To test this assumption 
here, we investigate the strain fluctuations using the non- 
affinity measured EE01 [25], T = J^((Su NA ) 2 ), where 
(5u NA = u u A denotes the nonaffine displacement of 
a cross-link and the brackets represent a network aver- 
age. This nonaffinity measure exhibits a cusp at the bend- 
ing rigidity percolation point, reflecting the criticality of 
the network's mechanics in this regime [6| [7], as shown in 
Fig. [2j). Furthermore, there appears to be a regime for 
sufficiently low k where r ~ L 2 independent k, lending 
credence to the basic assumption that Su lsA ~ L-f 12J. 




FIG. 3: (Color online) The shear modulus scaled with its 
affine prediction versus L scaled with Ana = £o/£t f° r vari- 
ous values of k in units of u£q. The open symbols indicate 
data ranges in the rigidity percolation regime where we ob- 
serve different scaling. The inset shows G scaled with |Az|^ 
as a function of \ Az\. 

Such nonaffine deformations imply a bending regime, 
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This prediction for the L-dependence of G is in fact born 
out by the numerical data, as shown in Fig. [2ji. This has 
the important implication that as L becomes large, Glc 
will exceed the affine shear modulus, and hence, an affine 
network deformation becomes more favorable. Thus, a 
NA-A transition fTUl IT4"l fT5] occurs when G LC ^ G A ; this 
is satisfied when L becomes comparable to the length scale 
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where = y/n/fi. Indeed, by plotting G/Ga as a func- 
tion of L/X NA we find a good collapse of the data to a 
universal curve, for which G/Ga ~ 1 when L/A NA > 1, 
as shown in Fig. [3] This strongly supports the existence 
of a NA-A transition driven by molecular weight in 3D 
fiber networks with connectivities well below Maxwell's 
CF isostatic point. In contrast, prior results for 2D net- 
works suggested A NA - i~ a , with a w 0.3 - 0.4 [T0"1HT] ; 
however, for such models it is unclear whether the NA-A 
transition is actually driven by fiber length, as for the 3D 
case presented here, or by the CF isostatic point [H] that 
coincides with the high-L limit for the Mikado model. 

nonlinear regime The length-controlled bending me- 
chanics also has important implications for the nonlin- 
ear elasticity of 3D fiber networks. Even in a bending 
dominated regime, stretching modes are excited at finite 
network deformations [T3], but to a higher order in the 
applied strain [7] [UJ [TB] . Specifically, assuming length- 
controlled nonaffine deformations, a transverse bend with 
an amplitude ~ 7L results in a stretch energy in the as- 
sociated bond, SEs ~ (J.e 2 £o, where e ~ (■yL/£ ) 2 + 0{^ i ). 
The onset of nonlinear network elasticity due to this ef- 
fect occurs at a strain 70, at which 8E$ becomes compa- 
rable to the leading order bending contribution, SEb ~ 
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FIG. 4: (Color online) a) The differential modulus scaled with 
its afSne prediction together with 70 and 7 A for networks with 
L — 18.3. b) Scaling predictions (solid lines) and numerical 
data for 70 (open symbols) and 7a (closed symbols). L is 
measured in units of to and k in units of pl\. The red and 
green points in the lower panel correspond to K = 10 -6 and 
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KL 2r y 2 /£ 3 . This stiffening saturates at a strain j A , set by 
the condition SEb + 6E$ ~ j?j 2 , at which the network's 
response becomes affine. Thus, the onset and completion 
of the stiffening regime are expected to scale as 
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These predictions are in excellent agreement with the nu- 
merical results over a broad range of L and shown 
in Fig. [4] Interestingly, the nonlinear threshold 70 corre- 
sponds to a stress <tq ~ k 3 ^ 2 L, which is distinct from the 
Euler buckling threshold [12] and the stiffening threshold 
in sub-isostatic spring networks [33] . 

The Phantom FCC model developed here, provides a 
powerful numerical model to probe the mechanics of 3D 
fiber networks with large system sizes. Using this model 
together with a scaling analysis, we have shown that even 
though the mechanical stability of 3D networks relies on 
the bending resistance of the constituent fibers, surpris- 
ingly for any k > 0, network mechanics becomes affine 
and independent of k when L > A NA [35]. This NA- 
A transition is induced because it becomes energetically 
unfavorable to accommodate length-controlled non-afhne 
deformations [12] at large molecular weights. Below this 
NA-A transition, in the length-controlled bending regime, 
we expect G ~ p 13 / 5 , for scmiflcxible polymers for which 
the cross-linking lengthscale is expected to scale with 
p -2 / 5 [21], or G ~ p 3 for stiff polymers, where p is the 
polymer length-density. These predictions may account 
for a recent report of G ~ p 2 ' es in collagen networks |17j. 

We conjecture that the results for the length-controlled 
regime also apply to models with additional interactions, 
other than fiber bending, that stabilize the network below 



the CF-threshold — including next-nearest neighbor inter- 
actions or bond-bending interactions for cross-links that 
fix a preferred bond-angle. Specifically, we expect that 
such networks exhibit an affine high-molecular weight 
limit even for arbitrarily weak additional interactions. 

This research was supported by the National Science 
Foundation under Grant No. NSF PHY05-51164 and by 
FOM/NWO. It's a pleasure to acknowledge discussions 
with E. Frey, M. Das, C. Heussinger, O. Stenull and T. 
C. Lubensky. We also thank OS and TCL for sharing 
their manuscript (Ref. |20| ) prior to publication. 



[9 

no 



[11 

[12; 

[13: 

[14 

[15: 

[ie: 
li- 
tis: 

[19 

[20 
[21 

[22: 

[23' 

[24 

[25' 



D. A. Fletcher and R. D. Mullins, Nature 463, 485 (2010). 
P. A. Janmey, S. Hvidt, J. Lamb, and T. P. Stossel, Nature 
345, 89 (1990). 

A. R. Bausch and K. Kroy, Nature Phys. 2, 231 (2006). 
P. A. Janmey et al., Nature Materials, 6: 48 (2007). 
J. C. Maxwell, Philos. Mag. 27, 294 (1864). 

C. P. B roedersz, X. Mao , T. C. Lubensky and F.C. MacK- 
intosh. |arXiv:1011.6535| (2011). 

M. Wyart, H. Liang, A. Kabla and L. Mahadevan, Phys. 

Rev. Lett. 101, 215501 (2008). 

A. J. Liu and S. R. Nagel, Nature 396, 21 (1998). 

M. F. Thorpe, J. Non-Cryst. Solids 57, 355 (1983). 

D. A. Head, A. J. Levine and F. C. MacKintosh, Phys. 
Rev. Lett. 91, 108102 (2003); Phys. Rev. E 68, 061907 
(2003). 

J. Wilhelm and E. Frey, Phys. Rev. Lett. 91, 108103 
(2003). 

C. Heussinger and E. Frey, Phys. Rev. Lett. 97, 105501 

(2006) ; C. Heussinger, B. Shaefer and E. Frey, Phys. Rev. 
E 76, 031906 (2007). 

P. R. Onck, T. Koeman, T. van Dillen and E. van der 
Giessen, Phys. Rev. Lett. 95, 178102 (2005). 
M. Das, F. C. MacKintosh and A. J. Levine, Phys. Rev. 
Lett. 99, 038101 (2007). 

A. R. Missel M. Bai, W. S. Klug and A. J. Levine, Phys. 
Rev. E 82, 041907 (2010). 

O. Lieleg et al., Phys. Rev. Lett. 99, 088102 (2007). 

S. B. Lindstrom, D. A. Vader, A. Kulachenko, and D. A. 

Weitz. Phys. Rev. E, 82, 051905 (2010). 

E. M. Huisman, T. van Dillen, P. R. Onck, and E. Van der 
Giessen, Phys. Rev. Lett. 99, 208103 (2007); E. M. Huis- 
man and T. C. Lubensky, Phys. Rev. Lett. 106, 088301 
(2011). 

G. A. Buxton and N. Clarke, Phys. Rev. Lett. 98, 238103 

(2007) . 

O. Stenull and T. C. Lubensky, (to appear). 

C. P. Broedersz and F. C. MacKintosh, Soft Mater 7, 3186 

(2011). 

B. A. Didonna and T. C. Lubensky, Phys. Rev. E 72, 
066619 (2005). 

M. Sheinman, C. P. Broedersz and F. C. MacKintosh (to 
appear). 

F. C. MacKintosh, J. Kas and P. Janmey Phys. Rev. Lett. 
75, 4425 (1995). 

For thermal semiflexible polymers we expect A NA ~ 



